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We report our theoretical results on the emergence of a partially-disordered state at zero temperature 
and its detailed nature in the periodic Anderson model on a triangular lattice at half filling. The partially- 
disordered state is characterized by coexistence of a collinear antiferromagnetic order on an unfrustrated 
honeycomb subnetwork and nonmagnetic state at the remaining sites. This state appears with opening 
a charge gap between a noncollinear antiferromagnetic metal and Kondo insulator while changing the 
hybridization and Coulomb repulsion. We also find a characteristic crossover in the low-energy excitation 
spectrum as a result of coexistence of magnetic order and nonmagnetic sites. The result demonstrates 
that the partially-disordered state is observed distinctly even in the absence of spin anisotropy, in marked 
contrast to the partial Kondo screening state found in the previous study for the Kondo lattice model. 

KEYWORDS: periodic Anderson model, geometrical frustration, triangular lattice, partial disorder, charge dis- 
proportionation, quantum critical point, partial Kondo screening 
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Rare-earth compounds have drawn considerable in- 
terest, since they exhibit fascinating behaviors, such as 
heavy quasi-particle formation, superconductivity, and 
non- Fermi liquid behavior, as a consequence of the in- 
terplay between localized electrons and conduction elec- 
trons.^^ The major outcomes of this interplay can be 
understood in terms of the two competing interactions, 
the Ruderman-Kittel-Kasuya-Yosida (RKKY) interac- 
tion and the Kondo coupling: The former favors magnetic 
ordering, whereas the latter promotes singlet for- 
mation between conduction and localized electrons.^' 
Their competition, in general, leads to a quantum criti- 
cal behavior between a magnetically ordered state and a 
Fermi liquid state, resulting in the rich phenomena in 
these systems. 

In the present study, we investigate the effect of geo- 
metrical frustration in the quantum critical region, with 
emphasis on the emergence of unusual orders. Experi- 
mentally, interesting coexistence of nonmagnetic Kondo 
singlet and magnetic order was observed in several rare- 
earth compounds with frustrated lattice structures, such 
as CePdAl^'^) and UNi4B.^°' Such coexistence was 
theoretically studied by the mean-field approximation 
(MFA) of a classical pseudomoment model. ^^'^^^ More 
recently, two of the authors and their collaborators 
showed a possibility of such coexistence, called the par- 
tial Kondo screening (PKS), by sophisticated variational 
Monte Carlo calculations for the Kondo lattice model 
(KLM) as well as the Kondo necklace model. 

While the existence of PKS is shown by explicitly 
taking into account quantum fluctuations and the spin- 
charge interplay, its origin has not been fully under- 
stood. In particular, it remains unclear whether the PKS 
appears in the absence of spin anisotropy. Moreover, the 
detailed nature of PKS, such as the electronic structure, 
was not clarified yet. The nature of phase transitions 
from the surrounding phases to the PKS state is also 
unclear. It is highly desirable to carry out the further 



analysis of the details of PKS phenomena. 

In this Letter, we investigate the ground state of a 
periodic Anderson model (PAM) on a triangular lattice. 
PAM is a parent model for KLM, and it includes charge 
degree of freedom in the localized level explicitly. The 
model allows a simple and straightforward MFA, which 
enables wide-ranging analysis compared to the previous 
numerical study. We will demonstrate that this itinerant 
electron model exhibits a partially-disordered (PD) state, 
a coexistence of a collinear antiferromagnetic (AF) order 
and nonmagnetic sites, without aid from spin anisotropy. 
The detailed nature of the PD phase and the associated 
phase transitions will be presented. 

We consider PAM whose Hamiltonian is given by^^^ 
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where c\^{cia-) and fl^{fia) are the creation (annihila- 
tion) operators of conduction and localized electrons at 
site i and spin a, respectively. The sum of (z, j) is taken 
over the nearest-neighbor sites on the triangular lattice. 
The first term in eq. (1) represents the kinetic energy of 
conduction electrons, the second term the onsite c-f hy- 
bridization between conduction and localized electrons, 
the third term the onsite Coulomb interaction for local- 
ized electrons, and the fourth term the atomic energy 
of localized electrons. PAM delivers KLM in the large 
U limit with one / electron per site. Hereafter, we take 
t = 1 as an energy unit and F > 0, and fix the electron 
density at half filling; n = ^^,^(4^^^^ + fLfio)/N = 2, 
where N is the total number of sites. We choose Eq = —1, 
while the following results are qualitatively the same for 
^0^0 within the noninteracting bandwidth. 

In order to obtain the ground state phase diagram 
of the model (1), we apply the standard Hartree-Fock 
approximation to decouple the Coulomb U term with 
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Fig. 1. (Color online), (a) Ground-state phase diagram of the 
model in eq. (1) at half-filling and Eq = —1. Schematic picture for 
the ordering pattern is shown for each phase. The size of the circles 
refiects the magnitude of local electron densities, and the arrows 
represent local spin moments. In the region for V < 0.2 and the 
narrow grey region between the F and AF metal, it is difficult to 
obtain a converged mean-field solution, (b) A three-site unit cell 
used in the Hartree-Fock calculation. Upper circles forming a tri- 
angular lattice and the dangling lower circles represent the c and 
/ electron degrees of freedom, respectively, (c) Schematic picture 
of the PD state. One of the three sublattices becomes nonmag- 
netic, while the remaining two sublattices forming a honeycomb 
subnetwork retain a collinear AF order. 
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Fig. 2. (Color online) V dependence of (a) energy gap, (b) mag- 
nitude of local spin moment and (c) local electron density at U = 2 
and Eq = —1. In (b) and (c), the data are plotted for c and / elec- 
trons at each sublattice separately, and the dashed curves represent 
the mean values of c and / contributions respectively. Note that 
site 2 and 3 are equivalent in each phase. 



preserving SU(2) symmetry in spin space. We employ 
three-site unit cell, as shown in Fig. 1(b), and calculate 
the mean fields by taking the sum over 200 x 200 grid 
points in the first Brillouin zone. The mean fields are de- 
termined self-consistently, by repeating calculations until 
they converge within the precision less than 10~^. 

Figure 1(a) shows the result for the ground-state phase 
diagram. The phase diagram is roughly divided into three 
regions; magnetic region for U > 2V^, nonmagnetic re- 
gion for [/ < 2V^, and intermediate region for /7 ^ 2V'^. 
Figure 2 shows the charge gap, magnitude of local spin 

mf(^) = yKFP + Kf)^+Kf^ is the ^ 

component of the spin operator for /(c) electron], and 

local electron density n{^''^ = {Y^a fLfi^i^a^i^)) ^ 
these three regions while varying V at /7 = 2. 

In the magnetic region for U > 2V^, a 120° AF order 
is dominantly stabilized. This state is a gapless metal, as 
shown in Fig. 2(a). The magnetic moments of conduction 
and localized electrons form 120° order in the entire re- 
gion of the AF metal, as schematically shown in Fig. 1(a) 
[see also Fig. 2(b)], except for a slight canting from the 
120° order seen for /7 > 6 and V < 0.4. The magnetic 
order is likely to be understood from the strong- KLM 
picture, in which the RKKY interaction should favor AF 
configuration between neighboring spins. In the smaller 
U region, the system becomes ferrimagnetic (F), whose 
magnetic state is schematically shown in Fig. 1(a), pre- 
sumably originating in the Stoner-like mechanism. 

On the other hand, in the regime for U < 21/^, the 
system becomes insulating [Fig. 2(a)] and nonmagnetic 
[Fig. 2(b)]. There, the c-f hybridization V gives rise to 
bonding and anti-bonding states, and the bonding states 
are fully occupied at half filling. This band insulator is 
called the Kondo insulator (KI). 

Between the magnetic and nonmagnetic regimes, the 



system exhibits a PD phase for U < 4.5. In this phase, 
one of the three sublattices becomes nonmagnetic, while 
the remaining two sublattices retain magnetic moments 
[Fig 2(b)]; see below for the detail. For larger [/, the 
system turns into an up-up-down type ferrimagnetic 
collinear state (uud) at ~ 5, and finally a 120° AF 
insulator for /7 > 5.5, as shown in Fig. l(a).^^^ 

Let us examine the nature of the intermediate PD 
state and the associated phase transitions. First of all, 
this state is insulating; the charge gap is nonzero, as 
plotted in Fig. 2(a). In the PD phase, the gap grows 
with increasing V, while its increase is suppressed for 
V > 0.98 after showing a shoulder-like behavior. We 
will return to this peculiar dependence later. Our result 
shows that the transition to the 120° AF metallic state is 
of first order with a jump of the gap, while it is continu- 
ous when entering to the Kondo insulator. We note that 
the results, especially the first-order transition between 
the noncollinear and collinear states, might be modified 
when taking account of a longer period order beyond the 
three-site unit cell, especially in the first-order transition 
between the noncollinear and collinear orderings. 

Second, the PD state exhibits peculiar coexistence of 
magnetic and nonmagnetic sites; the magnetic moments 
of both conduction and localized electrons vanish on one 
of three sublattices, while they remain finite on the re- 
maining two sublattices, as shown in Fig. 2(b). The mo- 
ment is larger for localized electrons than for conduc- 
tion ones, and the magnitude of each moment takes the 
same value between the two sublattices, as in the 120° 
AF metallic state for smaller V. The magnetic moments, 
however, show a collinear AF order on the two sublattice 
sites [Fig. 1(a)]. Therefore, the PD state is realized by 
relaxing the geometrical frustration; that is, the system 
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Fig. 3. (Color online). LDOS of conduction and localized elec- 
trons for each phase at [/ = 2 and Eq = —1. The data are taken at 
(a) V = 0.87, (b) 1.20, (c) 0.90, (d) 0.98, and (e) 1.10. (a) and (b) 
correspond to the 120° AF metallic phase and the KI phase, respec- 
tively, while (c)-(e) are for the PD phase. The thin vertical lines 
show the Fermi levels defined by the energy of the highest occupied 
level. Note that all the sites are equivalent in (a) and (b), and that 
the site 2 and 3 are equivalent in (c)-(e). {f)V dependence of the 
energy differences between the peak 3 and 4 and between the peak 
3 and 1. The dashed and dot-dashed curves denote 2Uml,^^ and 



2(3) 



"2(3)' respectively. For comparison, the energy gap in Fig. 2(a) 
is also plotted. 



becomes free from the magnetic frustration by sponta- 
neously self-organizing into the unfrustrated honeycomb 
AF network and the remaining nonmagnetic sites, as 
schematically shown in Fig. 1(c). Note that this peculiar 
state is obtained as the ground state, qualitatively dif- 
ferent from the entropically-driven PD states sometimes 
found in frustrated classical spin systems. 

Finally, the PD state is accompanied by charge dis- 
proportionation: the local charge density is higher at the 
nonmagnetic site than the AF sites, as shown in Fig. 2(c). 
Near the transition into the KI phase, the charge dispro- 
portionation decreases linearly with decreasing the dis- 
tance from the critical point, in contrast to the square- 
root behavior of the local spin moment in Fig. 2(b). This 
indicates that the collinear AF moment is a primary or- 
der parameter and the charge disproportionation is in- 
duced as a secondary effect. 

Now, we discuss the phase changes from the viewpoint 
of the electronic state. Figure 3 shows the local density of 
states (LDOS) for conduction and localized electrons for 
several V ditU = 2. In the AF phase, the energy spectrum 
is modulated due to the presence of static AF magnetic 
moments. As a result, there appear two sharp peaks, as 
shown in Fig. 3(a), whose energy separation is given by 
the mean- field Coulomb exchange energy 2Um{ . On the 



other hand, in the KI phase, there appear two peaks on 
both sides of the energy gap, as shown in Fig. 3(b). Their 
origin, however, is not magnetic but the c-f hybridization 
V. The peaks consist of the hybridized bands, which have 
a major contribution from the localized / levels. 

In the intermediate PD phase, the energy spectrum 
basically shows coexistence of the characteristic features 
of both the AF and KI phases [Figs. 3(c)-(e)]. At the 
nonmagnetic site (upper half in the figures), the LDOS 
shows two prominent peaks indicated by the peak 1 and 
2 in Fig. 3(c). They appear on both sides of the gap 
and are dominated by the contribution from localized 
electrons. This is characteristic of the hybridization ef- 
fect in Fig. 3(b). Meanwhile, at the magnetic sites (lower 
half), the LDOS exhibits two sharp peaks (peak 3 and 
4). The energy difference between the two peaks is close 
to 2Um^^^^^ as in the AF state, at least, in the region 
near the AF metallic phase, where the collinear AF mo- 
ment is well developed [see also Fig. 3(f)]. In addition, 
the energy difference between the peak 3 and 1 is scaled 
approximately to Um^^^y as plotted in Fig. 3(f). This 
indicates that the system gains energy to stabilize the 
PD state from the formation of the collinear AF order 
on the self-organized unfrustrated honeycomb network. 

With approaching the transition to the KI by increas- 
ing V, the gap-edge contributions at the nonmagnetic 
site (peak 1 and 2) smoothly vary to the peaks of the 
hybridized bands in the KI [Figs. 3(c)-(e)]. On the other 
hand, the LDOS at the AF sites changes its nature con- 
siderably with increasing V. As plotted in Fig. 3(f), the 
energy difference of the two peaks scales well to 2Um^^^^^ 
for small V as mentioned above, but it starts to deviate 
for y > 0.95: Finally, the two peak structures merge into 
the hybridized bands in the KI phase. This behavior is 
considered to be a crossover from the AF dominant to 
hybridization dominant regime. 

We further analyze the qualitative change of energy 
spectrum in the PD state in terms of the band struc- 
ture. Figure 4(a) shows the typical band structure in the 
PD phase. There are six bands (doubly degenerate each) 
under the three-sublattice ordering. The overall struc- 
ture does not change so much in the entire region of 
the PD phase, but an interesting qualitative change ap- 
pears near the Fermi level while increasing V . As shown 
in Figs. 4(b)- (d), although the highest occupied band 
(HOB) does not show remarkable change, the lowest 
unoccupied band (LUB) exhibits nontrivial change. For 
small V [Fig. 4(b)], the bottom of the LUB resides at 
the M' point in the Brillouin zone. This is attributed to 
the well-developed threefold AF order which opens a gap 
at the M' point. With increasing V, the energy at the F 
point gradually decreases relatively to the M' point, and 
becomes lower for V > 0.98 [Figs. 4(c) and 4(d)]. Thus, 
the gap structure changes from the indirect one between 
r and M' to the direct one at the F point. The latter 
smoothly turns into the hybridization gap in the KI.^^^ 
Hence, the result shows a crossover of the LUB nature 
from the AF dominant to hybridization dominant. 

The crossover of LUB is clearly reflected in the 
shoulder-like feature in the V dependence of energy gap 
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Fig. 4. (Color online). Energy dispersion at U = 2 and Eq = —1. 
The results are shown along the F-M'-K'-F line in the Brillouin 
zone, as shown in the inset of (a). In the inset, the solid hexagon 
represents the first Brillouin zone and the dashed one is for the 
folded Brillouin zone under three-sublattice ordering, (a) Overall 
structure of the energy dispersion at V = 0.98. (b)-(d) Energy dis- 
persion in the vicinity of the Fermi level at (b) V = 0.90, (c) 0.98, 
and (d) 1.10, respectively. The bold(blue) and bold-dashed(red) 
lines denote LUB and HOB. 



[Figs. 2(a) and 3(f)]. For V < 0.98 where the LUB is min- 
imum at the M' point, the indirect gap increases with V. 
Meanwhile, for V > 0.98 where the LUB bottom comes 
at the r point, the gap does not show substantial V 
dependence, as the energy difference between HOB and 
LUB is almost constant at the F point. 

The peculiar behavior of LUB implies rich possibilities 
caused by electron doping to the PD state. In particu- 
lar, in the crossover region, the LUB becomes extremely 
flat, as shown in Fig 4(c), which leads to a sharp peak of 
LDOS at the upper edge of the gap [Fig. 3(d)]. This sug- 
gests a possibility to have coexistence of heavy-fermion 
state and magnetic ordering, as experimentally observed 
in CePdAl.^^ It will be interesting to investigate such car- 
rier doping effect, possibly by more sophisticated method 
than the MFA. 

Finally, let us discuss our results in relation with the 
previous theoretical study for KLM.^"^^ As mentioned 
above, KLM corresponds to the large U limit of PAM 
with one / electron per site; the / electrons give local- 
ized moments, which couple with conduction electrons 
via the Kondo coupling J oc V"^ /U. Hence, the phase se- 
quence of AF, PD, and KI while changing V^/[/ in our 
results for PAM appears to well correspond to that of AF, 
PKS, and Kondo spin liquid (KSL) while changing J in 
KLM. The coincidence, however, needs careful consider- 
ation. The point is that our PD state takes place even in 
the absence of the spin anisotropy, whereas the PKS in 
KLM is hardly stabilized without spin anisotropy. Fur- 
thermore, our PD state appears in the relatively weak 
U regime, far from the strong U limit where KLM is 
justified. These results indicate that, besides the spin 
anisotropy, the charge degree of freedom via the c-f hy- 
bridization plays an important role in the stabilization 
of the PD state. 

When the spin anisotropy is introduced into the 



present PAM, the PD state is expected to extend to a 
wider range of parameters, and is presumably connected 
to the PKS state found in the strong-/7 limit, i.e., in 
KLM. Furthermore, we expect that the PD state in PAM 
is further stabilized by taking quantum fluctuations into 
account beyond the MFA, since the previous study indi- 
cated that quantum fluctuations stabilize PKS.^^^ Such 
extensions are left for future study. 

In summary, we have investigated the ground state of 
the periodic Anderson model at half filling on a triangu- 
lar lattice by the Hartree-Fock approximation. We found 
that the model exhibits a partial disorder with coexis- 
tence of magnetic order and nonmagnetic sites, which 
was not apparently observed in the related Kondo lat- 
tice model in the absence of spin anisotropy. We clarified 
the phase diagram in a wide parameter range and the na- 
ture of phase transitions from the surrounding phases to 
the partially-disordered phase. We also revealed the de- 
tailed nature of the partially-disordered state: it appears 
with gap opening and shows a characteristic crossover be- 
havior in the structure of the excitation spectrum from 
the dominantly-antiferromagnetic to dominantly-Kondo- 
insulating regime. 
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